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Here, we adapt the concept of transformational thermodynamics, whereby the flux of tempera- 
ture is controlled via anisotropic heterogeneous diffusivity, for the diffusion and transport of mass 
concentration. The n-dimensional, time-dependent, anisotropic heterogeneous Fick's equation is 
considered, which is a parabolic partial differential equation also applicable to heat diffusion, when 
convection occurs, for example in fluids. This theory is illustrated with finite element computations 
for a liposome particle surrounded by a cylindrical multilayered cloak in a water-based environment, 
and for a spherical multilayered cloak consisting of layers of fluid with an isotropic homogeneous dif- 
fusivity, deduced from an effective medium approach. Initial potential applications could be sought 
in bio-engineering. 
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INTRODUCTION 

This communication aims to target already existing 
systems which could enable the use of cloaking concepts 
in order to achieve control of three-dimensional processes 
using coated spheres consisting of concentric layers of ho- 
mogeneous isotropic diffusivity. Various applications al- 
ready implicate the use of concentric bilayered vesicles, 
one example being liposomes used for drug delivery 
Liposomes are concentric bilayered vesicles in which an 
aqueous volume containing a water-soluble drug is en- 
closed by a membranous lipid bilayer composed of nat- 
ural or synthetic phospholipids. One popular type of 
liposomes, known as the stealth liposomes Q are highly 
stable, long-circulating liposomes whereby polyethylene 
glycol (PEG) has been utilized as the polymeric steric 
stabilizer Q. Stealth and other liposomes use the con- 
cept of 'invisibility' in order to hide and evade the im- 
munosystem by coupling water-soluble polymers to the 
lipid heads. Therefore the polymer part of the molecule 
is dissolved in the aqueous environment, thus masking 
the liposomes from immune cells in the blood Other 
alternative applications to liposomes are nanoparticlcs 
based on solid lipids (SLN). These are composed of solid 
lipids stabilized with an emulsifying layer in an aqueous 
dispersion. This has benefits such as drug mobility. The 
release of the drug-enriched core of SLN is based upon 
Fick's first law of diffusion 

Another similar idea to the one presented in this paper 
is the concept of optical transparency resulting from the 
use of preparative reagents for the subcellular localization 
of fluorescently labelled tissues and organisms [f| . High 
resolution imaging techniques such as Laser Scanning Mi- 
croscopy (LSM) are able to provide high-resolution im- 
ages of biological samples. However the resolution of 
imaging whole organisms such as embryos by fluores- 



cently labelling certain components, becomes somewhat 
distorted due to the biological samples containing op- 
tically opaque regions. These opaque components are 
able to transmit, reflect, scatter as well as absorb light, 
which can lead to image distortions. Therefore certain 
commercially or non-commercially available reagents are 
available which are known as optical clearing reagents 
0, 0] and can render tissues and organisms transparent. 

Here we suggest a novel application to the fast-growing 
research area of cloaking, whereby a better control of 
light can be achieved through transformational optics, 
following the pioneering theoretical works of Pendry et al. 
[|| and Leonhardt to diffusion processes in biophysics. 
The aforementioned paper Q demonstrates the possibil- 
ity of designing a cloak that renders any object inside it 
invisible to electromagnetic radiation (using the covari- 
ant structure of Maxwell's equations), while the latter 
concentrates on the ray optics limit (using conformal 
mappings in the complex plane for Schrddinger's equa- 
tion). In both cases, the cloak consists of a meta-material 
whose physical properties (permittivity and permeabil- 
ity) are spatially varying and matrix valued. This route 
to invisibility is reminiscent of the work of Greenleaf et al. 
in the context of electrical impedance tomography fioj ]. 
Interestingly, the isomorphism between the anisotropic 
conductivity and thermostatic equations makes it pos- 
sible to control the pathway of heat flux in a station- 



ary setting, as observed by Fan et al. [ll| (see also [12 1 
for analogous cloaking in electrostatics) and experimen- 
tally validated by Narayano and Sato [l3[- However, time 
plays an essential role in diffusion processes, and manipu- 
lation of heat flux through anisotropic diffusivity requires 
greater care in a transient regime [l4l. Il5|. Interestingly, 
anisotropic diffusion is a well-known technique in com- 
puter vision [l(| aiming to reduce image noise without 
removing significant parts of the image content, typically 
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FIG. 1. Two-dimensional simulation for diffusion of chemi- 
cal species's concentration: Concentration is normalized to 1 
mol.m -3 on the left boundary with a flux boundary condition 
on right boundary with mass transfert coefficient of 5 m.s -1 , 
and symmetry boundary conditions on top and bottom; Two 
time point (t = 1 10 _6 s,(a,c) and t = 1.5 10 _5 s, (b,d)) simula- 
tions of mass diffusion in surrounding medium with diffusion 
constant of 2.1 10 _9 m.s -2 (C02-water) of a circular nano- 
size particle (nanobody) with a diameter of 1.5 10 _8 m of dif- 
fusion constant 1.9 10 _11 m.s" 2 (POPC-dehydrated). (b,d): 
Application of a cloak, surrounding the nanobody, which is 
of inner radius 1.5 10 _8 m and outer radius 3.0 10 _8 m and 
consists of 5 concentric layers. The first, third and fifth lay- 
ers from the inside of the cloak outwards have a diffusivity of 
4.586 10 _10 m.s -2 (sucrose in water) and respective thickness 
4.25 10~ 9 m, 5.25 l(T 9 m and 4 l(T 9 m. The second layer and 
the fourth layers have diffusitivities of 11 10~ 6 m.s~ 2 (Gaz 
ethanol-air) and 8.4 10" 9 m.s" 2 (liquid C02-Methanol) and 
identical thickness 7.5 10 _10 m. Note that in panels (b) and 
(d) isovalues of concentration (black curves) are bent around 
the nanobody, whilst they remain aligned outside the cloak. 



edges, lines or other details that may be important in the 



interpretation of the image 17]. Spatio-temporal differ 



ential equations of a reaction-diffusion type also appear 
in organogenesis models for the developments of limbs, 
lungs, kidneys and bones fl8j . 

The mathematical model described in this report is 
based upon Fick's laws of diffusion, derived by Adolf 
Fick in the year 1855 [l9|, which describes diffusion pro- 
cesses governing various contexts (conduction of elec- 
tricity, heat, concentration of chemical species etc, and 
even gray scale image in computer vision). Here, the 
diffusion coefficient is spatially varying (heterogeneous) 
and matrix valued (anisotropic). We show numerically 



that control can be achieved of three-dimensional pro- 
cesses (with a focus here on bio-engineering applica- 
tions), using coated spheres consisting of concentric lay- 
ers of homogeneous isotropic diffusivity, which mimick 
certain anisotropic heterogeneous diffusivity. Previous 
studies have only shown the control of diffusion processes 
with two-dimensional transformational thermodynamics 
13, llj]. Using a similar strategy, we extend earlier 
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works to three-dimensional diffusion processes and also 
discuss the issue of convection. 



MATERIALS AND METHODS 



Numerical simulations 



The time-dependent convection diffusion equation was 
implemented in the commercial finite element software 
COMSOL MULTIPHYS1CS. Two-dimensional computa- 
tions such as the one shown in Fig. [TJb) were performed 
using a laptop with 3Gb of RAM memory, but only re- 
quired 1Gb of RAM memory to achieve a complete con- 
vergence of the numerical solution to the parabolic prob- 
lem (the complete convergence was checked by refining 
the mesh size and the time steps, in the limit of the 3Gb 
of available RAM memory). For instance, the mesh for 
the fully converged numerical simulation shown in Fig. 
QJb) consisted of 16637 mesh points and 33142 trian- 
gular elements. This computation took less than two 
minutes using second order Lagrange equations and a 
time-dependent direct linear system solver (SPOOLES) 
with a symmetric matrix condition. Time step was 0.01s 
throughout the time interval [0s, 0.1s] with a relative tol- 
erance 0.01. However, three-dimensional computations 
shown in Fig. [21[3]were more computationally demanding 
than in the two-dimensional case. They were performed 
on a cluster of 16 computers with a total of 1Tb RAM 
memory. Nevertheless, fully converged numerical solu- 
tions required only 36 Gb of RAM memory to properly 
run (convergence of the results was double-checked refin- 
ing the mesh and time steps, what required up to 524Gb 
of RAM memory and 48 hours of computation) . For the 
results shown in Fig. [2][3l the mesh consisted of 152805 
mesh points and 905294 tetrahedral elements for the 3D 
computations, and computations took approximately 12 
hours (using up to 36Gb of RAM memory). More specif- 
ically, we used second order Lagrange elements and a 
congugate gradients time dependent linear system solver 
with an Algebraic multigrid (quality of multigrid hier- 
archy of 3) and time step of 0.01s throughout the time 
interval [0s, 0.1s], with a relative tolerance of 0.01. 
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RESULTS AND DISCUSSION 



Recent work used a change of coordinates in the 
time-dependent heat equation [20| to achieve a marked 
enhancement in the control of heat fluxes in two- 
dimensional media described by an anisotropic hetero- 
geneous conductivity (l4|. However, it has been known 
since the work of Fick 19[ that there is a deep anal- 
ogy between diffusion and conduction of heat or elec- 
tricity: Because of Fick's work, diffusion can be de- 
scribed according to the same mathematical formalism as 
Fourier's law for heat conduction, or Ohm's law for elec- 
tricity. Wc would like to use similar analogies between 
diffusion of heat and concentration of chemical species 
to propose an original strategy towards cloaking in bio- 
engineering/chemical engineering. A possible application 
is shown in Fig. [U through the creation of different lay- 
ers. This concept has already been used to a certain ex- 
tent in multivesicular liposomes which consist of bilayers 
of phosholipids. Using the diffusivity values published for 
a certain lipid-conjugated drug (1.9 10~ 11 m.s~ 2 for dehy- 
drated POPC) and coating these nanosized particles (30 
nanometers) with different layers of sucrose (diffusivity 
of 4.586 10" 10 m.s -2 ) and also thin layers containing sub- 
stances of higher diffusivities, results in the initial particle 
becoming invisible in the chemical environment (diffusiv- 
ity value of 2.1 10~ 9 m.s -2 ), which is used to represent 
an environment such as blood whose main constituent 



plasma is essentially composed of water [2l|, [22( . This can 
be seen with the comparison of the distribution of concen- 
tration without the cloak, Fig. DIa,c), an d with the cloak, 
Fig. [ljb,d). The consequences of the additional layers 
are more prominent at longer time points compared to 
shorter time points, t = 1.510~ 5 s (Fig. [TJ panels (c,d)) 
compared to t — 1 10~ 6 s (Fig. [U panels (a,b)). Thus, the 
effect of the different layers aids in maintaining high con- 
centrations of a substance in the center of the liposome 
for longer periods of time. This could have advantages 
in increasing drug stability for longer circulation times. 

Potential fabrication of this alternative liposome would 
involve additional procedures to the classical liposome 
fabrication steps currently used. More specifically, chlo- 
roform or chloroform-methanol or mixing with the lipid 
and hydrophobic organic solvents could still be used for 
the formation of vesicles. The appropriate removal of sol- 
vents involving rotary evaporation for extended time pe- 
riods could still be used. This could be followed by frozen 
storage before consequent steps. Secondly, replacing the 
classical hydration step by using an aqueous medium but 
with a higher concentration of sucrose compared to tra- 
ditional sucrose concentrations (with mixtures of glycine 
or alanine, which share the same diffusion coefficients 
and are also traditionally used for coating drugs), would 
allow the creation of sucrose layers during the forma- 
tion of micelles. Choosing appropriately sized vesicles 



would involve sonication techniques and analysis using 
techniques employed for vesicular structures such as dy- 
namic light scattering (DLS) equipment and transmission 
electron microscopy (TEM). 

Layers containing the gaseous phase or an appropri- 
ate replacement of a high diffusivity value (in the range 
of 10~ 6 to 10~ 5 m.s~ 2 ) would be the most complicated 
steps in the fabrication process. In practice, it would 
be more feasible to replace the gaseous layers used in the 
numerical stimulations by media of similar diffusivity val- 
ues, due to difficulties in initial fabrication and stability 
or maintenance of these layers. 

It should be noted that the examples of specific sub- 
stances chosen for the simulations can easily be replaced 
by alternative appropriate substances with the same dif- 
fusion coefficient values. In addition, if potential cloaking 
applications are non-objectionable to the use of chloro- 
form, it should be noted that this replacement has al- 
ready been calculated to be a good substitute for the su- 
crose layers (layer one, three and five). This is shown in 
Supplementary Fig. 1, wherein there is improved cloak- 
ing as the isovalue curves for concentration outside the 
multilaycred structure are nearly perfectly aligned (pan- 
els (b),(d)) in contradistinction to what can be observed 
in panels (a), (c) for a nanoparticle not surrounded by 
a cloak. Supplementary Fig. 2 shows that the max- 
imum concentration within the nanoparticle is always 
lower when it is surrounded by the cloak (and that its 
variation is dramatically reduced because of the cloak). 
This is demonstrated by calculating the concentration 
at all points along a line passing through the center of 
the nanoparticle without (Supplementary Fig. 2, panels 
(a),(c)) and with (Supplementary Fig. 2, panels (b),(d)) 
a surrounding cloak. Note that this is achieved with sim- 
ply 5 concentric layers, 3 of which have same diffusivity. 
A similar type of profile for concentration can be observed 
in Supplementary Fig. 4, for a three-dimensional cloak 
with 20 layers described in Supplementary Fig. 3. 



n-dimensional transformed convection-diffusion 
equation 

We consider the convection-diffusion equation which is 
a parabolic partial differential equation combining the 
diffusion equation and the advection equation. This 
equation describes physical phenomena where particles 
or energy (or other physical quantities) are transferred 
inside a physical system due to two processes: diffusion, 
which results in mixing and transport of chemical species 
without requiring bulk motion (it is a random walk of 
particles/molecules towards certain equilibrium state i.e. 
homogeneous distribution of chemical species inside a re- 
gion), and convection, whereby collective movements of 
ensembles of molecules take place (usually in fluid) which 
in essence use bulk motion to move particles from one 
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place to another place [23| . In its simplest form (when the 
diffusion coefficient and the convection velocity are con- 
stant and there are no sources or sinks), the convection- 
diffusion equation in a domain f2 (with a chemical source 
outside) can be expressed as [24j 



dc 
dt 



E 



E 



_d_ 



■ViC , 



(1) 



where c represents the mass concentration (in biochem- 
istry) evolving with time t > 0, k is the chemical diffu- 
sion in units of m 3 .s _1 . and v the velocity field. We note 
that Fick's equation is written in a general form, where 
x = (x±, ...,x n ) is a variable in an n-dimensional space. 
Accordingly, sums stretch from i,j = 1, ...n (here appli- 
cations are sought in 2d and 3d spaces, so n = 2 or 3). 
It is customary to put matrix k in front of the spatial 
derivatives when the medium is homogeneous. However, 
here we consider a heterogeneous (possibly anisotropic) 
medium, hence the spatial derivatives of k might suf- 
fer some discontinuity (mathematicall y, p artial deriva- 
tives are taken in distributional sense [25j . hence trans- 
mission conditions ensuring continuity of the heat flux 
kVc are encompassed in (p}). Physically, the diffusion 
flux — kVc measures the amount of substance that will 
flow through a small volume during a short time interval 
(mol.m~ 3 .s^ 1 ). 

Upon a change of variable x = (xi, X2, X3) — > y = 
(2/1,2/2,2/3) described by a Jacobian matrix J such that 
Jij = dyi I dxj , (JTJ) takes the form: 



detJij dt 



' C = E^( 



1 



dc 

J ij C 3 



5 detJij ' 



(2) 



where k[_ = jKj T det(J) 1 and vL = det(J) 1 3 T v_ are the 
transformed diffusivity and velocity, respectively. 



Jacobian matrix and transformed diffusivity for 
cloaking 

Let us now consider the following transform [Io| 

F(x)=(l + hx\]-^-, (3) 



where | x \= \J x\ + x\ + . . + x 2 . This function is 
smooth except at point O = (0, ..,0). It blows up the 
point O to the hyperspherc of radius | x |= I, while 
mapping the hypersphere of radius | x \= 2 to itself. 
Moreover, F(x) = x at the boundary | x |= 2. 
Defining the Jacobian matrix J as Jij = dFi/dxj, we 
find: 



1 + 5 I « 



I 



1 



(4) 



where I is the n x n identity matrix and x = xj \ x |. 
This Jacobian is well defined everywhere except at x = 0. 
We note that J is symmetric, x is an eigenvector with 
eigenvalue 1/2 and x J ~ is a n — 1 dimensional eigenspace 
with eigenvalue l/2 + l/|x|ina space of dimension n. 
The determinant of the Jacobian follows easily: 



det(J) = Ul + \ I x 



(I x I +2)"- 
2 n I x I 11 " 1 



(5) 



There are two cases of practical importance: n = 2, for 
which writing r =| x \= \J x\ + x\ = 2(| y \ —1) = 



2(r' — 1), one can see that the eigenvalues of the matrix of 
transformed diffusivity k' behave like r and r _1 as r — > 0; 



3 for which writing 



x\ 



■I' 



n 

2(1 y I ~ 1) = 2(r' — 1), one can see that the matrix 
of transformed diffusivity k' has one eigenvalue which 
behaves like r 2 and two like r° as r — > 0. Interestingly, 



for n > 4, one eigenvalue of k' behaves like 



and the 



remaining n — 1 behave like r n ~ . This shows that only 
the case n = 2 leads to a singular matrix of diffusivity k' 
at the inner boundary of the cloak (the circumferential 
eigenvalue becomes infinite) , a fact already noticed in the 
context of cloaking for electric impedance tomography 
[Tol 26 1 . However, this matrix is always degenerate at 



the inner boundary of the cloak, irrespective of the space 
dimension. Similarly, j/ is a null vector on the inner 
boundary of the cloak only in space dimension 2. This 
analysis provides evidence that spherical cloaks should 
be easier to construct than circular cloaks. 
The parameters of the bio-cloak need to be further anal- 
ysed in polar and spherical coordinates in order to sim- 
plify the cloak's design and the numerical implementa- 
tion. 



On the choice of reduced parameters for a bio-cloak 
in polar and spherical coordinates 

We first note that if we multiply both sides of ^ by 
detJij and let detJij inside the partial space derivatives, 
we retrieve the usual form of the convection-diffusion 
equation, albeit with anisotropic coefficients. We realise 
this is not legitimate in general as by doing so we add 
an extra term £\ ^. JijKij(y) J^^-(detJij^-) in ©, but 
we numerically checked that this term can be sufficiently 
small that it does not significantly affect the solution of 
the original transformed equation (J2J). Physically, this 
manipulation results in preserving the direction of the 
diffusion flux —k/Vc (since detJij is a scalar). However, 
it affects its continuity (since detJij is heterogeneous). 
Such a manipulation is known in the transformational 
optics community to lead to transformed equation with 
reduced parameters [27}. We now observe that from the 
function F(r) = R\ + r(i?2 —Rxj/Ri counterpart of |3]), 
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wherein R\ = 1 and i?2 = 2, in polar (r,8) (resp. spheri- 
cal (r,6,(j>)) coordinates, which blows up a point O to the 
disc (resp. the sphere) of radius i?i and maps the disc 
(resp. the sphere) of radius i?2 to itself in polar (r',0') 
(resp. spherical (r', 6>', </>')) coordinates, the transformed 
diffusivity can be expressed for a cylindrical cloak as: 



i? 2 \Vr'-i?i x2 



i?2 — Rl , 
and for a spherical cloak as 



, Kg, 



B 2 



i?2 — -Rl 



(6) 



R 2 Yfr'-R^ 4 



R2-R1J \ r 

R 2 Vfr'-R^ 2 

Ka, — K j,, 



i?2 — -Ri 



(7) 



thermal cloaks in 

Kp, 



where R\ and R2 are the interior and the exterior radii of 
the cloak. One note that when r' tends to R\ , k' t , goes to 
zero and Kg, remains constant in ([6]), whereas n' r ,, Kg, and 
k'^i all go to zero in ([7|). This means that thanks to the 
reduced coefficients the matrix of transformed diffusivity 
k' now has one eigenvalue which behaves like r 2 and one 
like r° as r — > in the cylindrical case i.e. we no longer 
have an eigenvalue which blows up on the inner boundary 
of the cloak. Likewise, we now have one eigenvalue which 
behaves like r 4 and two like r 2 , instead of one behaving 
like r 2 and two like r° when r tends to zero, if we use 
reduced parameters in the spherical case. Thus, reduced 
parameters are an obvious choice in the cylindrical case, 
since Kg, is a constant in ([6]), and were implemented for 

, [15J. However, in spherical case, 

are no longer constant in ([7]). Nevertheless, 
we need use such reduced coefficients to get rid of the 
coefficient sitting in front of the time derivative in ([2]). 
These heterogeneous anisotropic (reduced or not) pa- 
rameters can be approximated by piecewise constant 
isotropic coefficients, making use of an effective medium 
approach, as detailed in the supplementary material, 
which justifies the implementation of multilayered cylin- 
drical and spherical cloaks with concentric isotropic ho- 
mogeneous thin layers. The choice of reduced parameters 
led in |14| to a multilayered thermal cloak with piece- 
wise constant and high contrast diffusivity, and we refer 
to values of diffusivity and computations therein for the 
cloaking effect for a concentration of chemical species in 
2D with 20 layers (there is a one-to-one correspondence 
between Fourier's heat equation and Fick's equation, the 
unknown being either the temperature or the concentra- 
tion). However, Fig. [1] of the present paper clearly shows 
cloaking can be achieved with simply 5 layers with mod- 
erate contrast in diffusivity. Moreover, in the spherical 
case, the choice of reduced parameters leads to a differ- 
ent set of parameters (see supplemental material) for the 
multilayered bio-cloak, which we now study numerically. 



Numerical illustation 



For illustrative purposes, we focus here on a spherical 
cloak consisting of 20 concentric layers with diffusitivity 



ranging from 2.5 10 m 



m .s 



further 



details can be found in the supplemental material. In 
order to emphasize the power of the approach, we now 
consider a cloak substantially larger than that in Fig. [1] 
Indeed, as observed in [!§] in the context of the trans- 
formed heat equation, if we rescale the coordinate system 
as x — t- sx and the time variable t — > s 2 t with a dimen- 
sionlcss factor s, the transformed Fick's equation © re- 
mains unchanged, provided we assume velocity is ruled 
out. More precisely, the cloak of Fig. [T]has been scaled 
up by a factor s = 10 2 i.e. its inner radius 1.5 10 _6 m and 
its outer radius 3.0 10" 6 m in Fig. [2] Accordingly, time 
should be scaled up by a factor s 2 = 10 4 . We show the 
distribution of concentration for a selectin of time points 
ranging from t = 5 10" 3 s to t = 2.5 10~ 2 s in Fig. O 
Moreover, in Fig. [2j the cloak is in presence of a chemi- 
cal specy with concentration normalised to 1 mol.m for 
simplicity (taking any other concentration C will simply 
lead to a color scale in Fig. [2] and [3] ranging from to 
C mol.m -3 ), which is set on the right-hand side of the 
computational domain (a cube of sidelength 8.0 10 _6 m). 
On the opposite (left-hand) side, we set the usual flux 
condition (kVc) • n = iVo + fc c (c;, — c) (with c the, as yet 
unknown, solution to (fl]) and n the unit outward normal 
to each side of the cube), where k varies within the range 
given in the caption of Fig. [2] inside the layers of the 
cloak (see supplemental material for more details), and 
k = 1 m 2 .s _1 in the inner core and outside the spherical 
cloak and k c is the mass transfer coefficient, iVo is the 
inward flux and c\, is the bulk concentration of chemical 
species in the cubic domain. Here, we consider k c = 5 
mol.s -1 , No — mol.m -2 . s _1 and q, = mol.m" 3 . Fi- 
nally, we set insulation (or equivalent symmetry) condi- 
tions (kVc) • n = on the four remaining sides of the 
cube. More details on implementation of such a diffu- 
sion model in finite elements may be found in (28j . This 
phenomenological model of a bio-cloak exhibits the fol- 
lowing features: The concentration of chemical species 
nearly vanishes inside the inner sphere of the bio-cloak 
at time t = 0.005s, see Fig. [DJa). In the optical setting, 
such a sphere is called invisibility region, as no scattering 
obstacle placed inside this region can be detected 
In a bio-physical setting, this zone acts as a protection 
from any potential chemical attack. However, the con- 
centration of chemical species increases steadily over time 
until it reaches half the value of the concentration which 
is set on one side of the cubic computational domain, see 
Fig. EJJd). We note that the concentration is always uni- 
form inside this inner sphere at any time point, see Fig. 
[5Ja)-(d). Such a bio-cloak therefore offers some kind of 
protection from attack of chemical species since the con- 
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ccntration is uniform in its inner sphere at any time point 
and the concentration therein at any time point is always 
smaller than it would be without a cloak (see supplemen- 
tal material). Moreover, it never exceeds a value, which 
in our configuration is half the applied concentration (it 
can be seen that this is due to the fact that the cloak has 
its origin in the center of the cubic domain). 



(a) (b) 




FIG. 2. 3D plot of concentration [mol/m 3 ]: (a) t=0.005s; 
(b) t=0.01s; (c) t=0.015s; (d) t=0.025s. It has been checked 
that 3D plots are as in (d) for t > 0.025s (steady state). 
Spherical cloak of inner radius 1.5 10 -6 m and outer radius 
3.0 10 _6 m consists of 20 concentric layers with diffusity rang- 
ing from 2.5 10~ 6 m 2 .s~ 1 to 1.7 10 _2 m 2 .s _1 . The core and 
outer medium have same diffusivity 1.5 10 _5 m 2 .s _1 . 



(a) (b) 




FIG. 3. 2D plot of concentration [mol/m 3 ] corresponding to 
a slice of 3D plot in Fig. [2] in the horizontal plane passing 
through the center of the cloak: (a) t=0.005s; (b) t=0.01s; 
(c) t=0.015s; (d) t=0.025s. It has been checked that 2D plots 
are as in (d) for time steps t > 0.025s (steady state). 



CONCLUSION 

In conclusion, in this brief communication we introduce 
a coordinate transformation approach to control diffusion 
processes via anisotropy with an emphasis on concen- 
tration of chemical species for potential applications in 
bio-physics or bio-enginccring. Not only does the form 
of the transformed convection-diffusion equation involve 
an anisotropic heterogeneous diffusivity, but it also re- 
quires a spatially varying coefficient in front of the time 
derivative, as well as an anisotropic heterogeneous ve- 
locity held. In order to be able to design a structured 
cloak for such diffusion processes, we have simplified this 
transformed equation by introducing so-called reduced 
coefficients, which preserve the direction of the diffusion 
flux (but create an impedance mismatch between the 
cloak and surrounding medium) and by further assum- 
ing a small velocity in a homogenization approach. Our 
theoretical results are exemplified with numerical simula- 
tions in two- and three-dimensional settings. It should be 
emphasized that the issue of convection, which is of par- 
ticular importance for diffusion in fluids, such as in living 
organisms, will require further comprehensive studies. 
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tion de La Recherche Medicale. 



[1] Martins, S., Sarmento, B., Ferreira, D. C, Souto, E. B. 
2007 Lipid-based colloidal carriers for peptide and pro- 
tein delivery-liposomes versus lipid nanoparticles. Int J 
Nanomedicine. 2(4): 595-607. 

[2] The Danish National Research Foundation 
for Biomembrane Physics (MEMPHYS) at 
the University of Southern Denmark, Odense. 
Bioimaging, Science in Your Eyes. 2008 
http://www.scienceinyoureyes. com/index. php?id=80 

[3] Barrajon-Catalan, E., Menendez-Gutierrez, M.P., Falco, 
A., Saceda, M., Catania, A. and Micol, V. 2011. Immuno- 
liposomes: A Multipurpose Strategy in Breast Cancer 
Targeted Therapy in Breast Cancer - Current and Al- 
ternative Therapeutic Modalities (Ed. E. Gunduz and M. 
Gunduz), ISBN 978-953-307-776-5, InTech. 

[4] Elbayoumi, T.A., Torchilin, V.P. 2006 Enhanced accu- 
mulation of long-circulating liposomes modified with the 
nucleosome-specific monoclonal antibody 2C5 in various 
tumours in mice: gamma-imaging studies. Eur. J. Nucl. 
Med. Mol. Imaging 33(10):1196-1205. 

[5] Muller, R.H., Dingier, A., Weyhers, H. Muller, R.H., zur 
Muhlen, A. 1997 Feste lipid nanopartikel (SLN). In Phar- 
mazeutische Technologie: Moderne Arzneiformen, (Ed. 
R.H. Muller, G.E. Hildebrand), pp. 265-272. Stuttgart: 
Wissenschaftliche Verlagsgesellschaft 

[6] Hama H., Kurokawa H., Kawano H., Ando R., Shimogori 
T., Noda, H., Fukami, K., Sakaue-Sawano A., Miyawaki 
A. 2011 Scale: a chemical approach for fluorescence imag- 
ing and reconstruction of transparent mouse brain. Nat. 



7 



Neurosci. 14(11):1481-1488 [20 
[7] Seisenberger, G., Ried M.U., Endress, T., Buning H., 

Hallek, M., Brauchle, C. 2001 Real-Time Single-Molecule [21 
Imaging of the Infection Pathway of an Adeno- Associated 
Virus, Science 294 (5548): 1929-1932 [22 
Pendry, J.B., Schurig, D. and Smith, D.R. 2006 Control- 
ling Electromagnetic Fields, Science 312, 1780 
Leonhardt, U. 2006 Optical Conformal Mapping, Science 
312, 1777 

Greenleaf, A., Lassas, M. and Uhlmann, G. 2003 On 
nonuniqueness for Calderon's inverse problem, Math. [23 
Res. Lett. 10, 685-693 

Fan, C.Z., Gao, Y. and Huang, J. P. 2008 Shaped graded [24 
materials with an apparent negative thermal conductiv- 
ity Appl. Phys. Lett. 92, 251907. [25 
Nicorovici, N.A.P., Milton, G.W., McPhedran, R.C. 
and Botten, L.C. 2007 Quasistatic cloaking of two- 
dimensional polarizable discrete systems by anomalous [26 
resonance, Opt. Express 15, 6314-6323. 
Narayana, S. and Sato, Y. 2012 Heat flux manipulation 
with engineered thermal materials, Phys. Rev. Lett. 108, [27 
214303. 

Guenneau, S., Amra, C. and Veynante, D. 2012, Trans- 
formation thermodynamics : cloak and concentrator for 
heat, Opt. Express 20(7): 8207-8218. [28 
Schittny, R., Kadic, M., Guenneau, S. and Wegener, M. 
2012 Experiments on transformation thermodynamics: [29 
Molding the flow of heat l|arXiv:1210.2810|l 
Hallinan, P.L., Gordon, G.G., Yuille, A.L., Giblin, P. and 
Mumford, D. 1999 Two and three dimensional patterns 
of the face, Natick, Massachusetts: A.K.Peters. [30 
Perona, P. and Malik, J. 1990, Scale-space and edge de- 
tection using anisotropic diffusion, IEEE Transactions on 
Pattern Analysis and Machine Intelligence 12(7): 629- [31 
639. 

Kondo, S. and Miur, T. 2010 Reaction-diffusion model 
as a framework for understanding biological pattern for- 
mation. Science 329: 1616-1620. 

A. Fick 1855 On Liquid Diffusion. Philosophical Maga- 
zine 4(10): 30-39. 



Fourier, J. 1822 Theorie analytique de la Chaleur, Paris: 
Firmin-Didot pere et fils. 

Martinez, Isidoro. Mass Diffusivity D ata. 1995-2013. 
|http: //webserver.dmt.upm.es/ ~isidoro/| 
Daykin, C.A. , Savage, A. K., Wulfert, F. 2013 
Univeristy of Nottingham and Engineering and 
Physical Sciences Research Council (EPSCR). i- 
metabolomics: spectra and diffusion coefficients 
of common metabolites found in blood plasma. 



http://www.nottinghamcourses.com/pharmacy/documents/i-metabo 
Stanford, A.L. 1975 Foundations of biophysics, New- 
York, London: Academic Press 

Cussler, E.L. 1997 Diffusion Mass Transfer in Fluid Sys- 
tems, New- York: Cambridge University Press. 
Jikov, V.V., Kozlov, S.M. and Oleinik, OA., 1994 Ho- 
mogenization of Differential Operators and Integral Func- 
tional, New- York: Springer- Verlag 

Kohn, R.V., Shen, H., Vogelius, M.S. and Weinstein, 
M.I. 2008 Cloaking via change of variables in electric 
impedance tomography, Inverse Problems 24: 015016. 
Schurig, D., Mock, J. J., Justice, B.J., Cummer, S.A., 
Pendry, J.B., Starr, A.F. and Smith, D.R. 2006 Meta- 
material electromagnetic cloak at microwave frequencies. 
Science 314: 977-80. 

Morton, K.W. 1996 Numerical Solution of Convection- 
Diffusion problems, London: Chapman and Hall 
Wolfe, C.A, James, P. A, Mackie, A, R., Ladha, S., Jones, 
R. 1998 Regionalized Lipid Diffusion in the Plasma Mem- 
brane of Mammalian Spermatozoa. Biology of reproduc- 
tion 59, 1506-1514. 

Gaede, H.C and Gawrisch, K. 2003 Lateral Diffusion 
Rates of Lipid, Water, and a Hydrophobic Drug in a 
Multilamellar Liposome. Biophys J. 85 (3): 1734-1740. 
Kwong P.D., Wyatt R., Robinson J., Sweet R.W., So- 
droski J., Hendrickson, W.A. 1998 Structure of an HIV 
gpl20 envelope glycoprotein in complex with the CD4 
receptor and a neutralizing human antibody. Nature 393 
(6686): 648-659 



